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Abstract 

Qj • Based on the law of mass action (and its microscopic foundation) and mass 

1^ ■ conservation, we present here a method to derive consistent dynamic models for 

the time evolution of systems with an arbitrary number of species. Equations 
are derived through a mechanistic description, ensuring that all parameters have 

^ I ecological meaning. After discussing the biological mechanisms associated to the 

logistic and Lotka-Volterra equations, we show how to derive general models for 



O 



O I trophic chains, including the effects of internal states at fast time scales. We show 

Q ! that conformity with the mass action law leads to different functional forms for 






the Lotka-Volterra and trophic chain models. We use mass conservation to recover 



^>j, the concept of carrying capacity for an arbitrary food chain. 
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1. Introduction 

There exists a multitude of models for trophic interactions. For example, 
Royama (1971) and May (1974) describe different alternatives to model the same 
interactions, and Berryman et al. (1995a) give a table of twenty five alternatives 
to model predator-prey systems. However, the calibration and validation of these 
models with experimental and observational data is systematicaly lacking, and 
most ecologists prefer to adjust time series data with empirical models that have 
no connection to the specific ecological processes (Solow, 1995). 

For single-species population dynamics, the logistic equation is the basic 
paradigm, introduced in almost any ecology textbook. It accurately predicts pop- 
ulation densities in systems such as bacterial batch cultures (Schlegel, 1992) and 
human populations (Banks, 1994) and, when generalized, describes the dynamics 
of many single species populations in both laboratory and field (Gause, 1934; AUee 
et al., 1949; Thomas et al., 1980; Berryman and Millstein, 1990). It is applica- 
ble to multiple situations in ecology and biology (Banks, 1994) and bioeconomics 
(Clark, 1990). However, the logistic equation has been criticized with the argu- 
ment that the underlying carrying capacity concept has no mechanistic meaning, 
being simply a fitting parameter (Kooi et al., 1998), and obscuring the relation 
between population growth and resource availability (Getz, 1984). 

For multi-species population dynamics, the basic model is the Lotka-Volterra 
equations. They are the basis of almost all the theory of trophic interactions. A 
further development of these equations was the recognition that there exist limits 
to the capacity for consumption, leading to the introduction of the HoUing func- 
tional response, commonly called Type II (HoUing, 1959). It was later verified that 
the functional form of this curve coincides with the Monod function used in mi- 
crobiology and the Michaelis-Menten mechanism in enzyme kinetics. The HoUing 
Type II functional response lies at the heart of current trophic chain dynamics 
theory (Oksanen et al., 1981). 

In 1928 Volterra adopted the mass action principle of chemical kinetics to 
write the dynamic equations for the densities of a prey-predator system (Berryman, 
1992). In 1977, Nicolis and Prigogine showed that the logistic equation could be 
derived in analogy with chemical kinetics, using the mass action law and a mass 
conservation principle. The mass action law lies at the heart of most population 
dynamics theory, as in epidemiology (Anderson and May, 1991) and in structured 
population models (Metz and Diekmann, 1986; Metz and de Roos, 1992). 

Mass conservation is a controversial issue in population dynamics. Some au- 
thors have argued that population dynamics models do not have to conform to 



mass conservation (e.g. Berryman et a/., 1995b). However, the effect of mass 
conservation on the dynamics of communities has been ascertained (De Angehs et 
a/., 1989). 

On the other hand, one of the problems with the complexity of ecological 
systems is that there may exist internal states of the systems which we cannot 
measure (Arditi and Ginzburg, 1989). There is a systematic way of eliminating 
these variables from the system description, if their dynamics occur at faster time 
scales than the time scale of population dynamics. In physics this is called the 
adiabatic approximation (Haken, 1983), and in chemistry it is called the quasi- 
steady state assumption (Segel, 1988; Segel and Slemrod, 1989; Borghans et ai, 
1996; Stiefenhofer, 1998). This approach has been used in ecology to distinguish 
between different time scales (O'Neill et a/., 1986; Michalski et al., 1997). 

In this paper, we take the chemical kinetics analogy to its full consequences, 
showing how to derive the population dynamic equations of an arbitrary food web 
from the ecological mechanisms of interaction. This approach is based on the fact 
that both organisms and molecules are discrete entities that interact with each 
other. The advantage of this analogy is that we can bring from physics to ecol- 
ogy the knowledge of statistical mechanics about the transition from individual 
motions to macroscopic behavior, deriving the precise limits of validity of the de- 
terministic population dynamics description (Maurer, 1998). This program unifies 
the deterministic with the stochastic approach to population dynamics, as far as 
local densities of individuals do not fluctuate too much around the average density 
of the whole population. 

This paper is organized as follows. In the next section we review the main 
techniques for deriving evolution equations of chemical kinetics, with the necessary 
modifications for ecological systems. The general evolution equations are taken 
in accordance to the mass action law. We then impose a conservation law which 
is equivalent to the assumption of closedness of the ecological system. This mass 
conservation law makes it possible to model a renewable resource, leading to the 
concept of carrying capacity. To account for internal states occurring at fast time 
scales, we introduce the mechanism that leads to a Michaelis-Menten resource 
uptake, which can be compared with logistic type mechanisms. In section 5, we 
show that the correct application of the mass action law leads to a new Lotka- 
Volterra type system of equations, for two-species and n-species interactions. In 
section 6, we derive the general form of the evolution equations for a trophic chain 
where species can have internal states. Taking limits of ecologically significant 
parameters, more complicated models are reduced exactly to simpler ones. 



2. From chemical kinetics to ecological mechanisms: the mass ac- 
tion law 

In order to construct population dynamics models with clear and ecologically 
significant mechanisms and precise conditions of validity and applicability, we now 
introduce the analytical technique that will be used in this paper. 

At the scale of interatomic distances, the motion of molecules in a solution 
is random. When two eventually binding chemical species collide, a new molecule 
appears, decreasing the mole number of the initial chemical species and increasing 
the mole number of the newly formed chemical species. Analogously, in population 
dynamics, if we assume randomness in the motion of individuals, the interaction of 
individuals with a resource is a collision. At collision, the individual can consume 
the resource — binding — or simply ignore it. Therefore, both systems can be 
considered similar and, at the macroscopic level, the mean densities (mole number 
per unit volume, in chemical kinetics; number of individuals per unit area or 
volume, in population dynamics) are described by the same evolution laws. 

So, we consider a closed area (or volume) 5" with several species or resources, 
Aj, j = 1, . . . , m, with number of individuals given by rij. Interactions in S are 
described by the collision diagrams 

Z/iiAi H hUimAra^''' fJ^ilAi^ h /U^^A^ , z = 1, . . . , n (2.1) 

where Uij and Hij are positive parameters measuring the number of individuals 
that are consumed and produced, with Uij being integers, and the constants r^ 
measure the rate of the interaction. 

Suppose further that the species Aj is well distributed in S with mean density 
ttj = Uj/S. It follows from a master equation approach (van Kampen, 1992, pp. 
166-172) that 

^ = E^^(^^.-^..)«r---C'", j = i,---,m (2.2) 

1=1 

Equation (2.2) expresses the mass action law and is derived using the following 
assumptions, establishing its limits of validity, 
i) At each instant of time, the densities of each species are approximately con- 
stant over the finite territory S. 
ii) The densities are low. 

iii) Individual motions are independent of each other, in such a way that the 
collision frequency is proportional to the product of the probability densities 
of finding the different individuals in a small region. 



iv) Interaction probabilities are independent of the past history of organisms. 
v) The motion of individuals is random ans is due to some form of collision with 

the environment, as in Brownian motion (Haken, 1983). 

Adopting this analogy between ecological and chemical interactions, it is pos- 
sible to derive the usual interaction laws found in the ecological literature, with 
the advantage that now these evolution laws have a precise mechanistic meaning 
given by the collision diagrams (2.1). 

In the following, we will make an additional simplification. Defining the order 
of a reaction as Yl7=i ^ij^ ^^ *^^^ only consider second or lower order reactions, 
^•^•' Yl7=i ^ij — 2' since higher order reactions have comparatively negligible prob- 
abilities of occurrence. 

In the following and to simplify the notation, we will represent species and 
species densities in diagram (2.1) and equations (2.2) by the same symbol Aj. 

3. Logistic autotrophs 

Let us represent by A^ the density of individuals of a species per unit of area 
or volume. Let us represent resources by A. Schematically we can represent the 
ecological interaction — species consuming resources and reproducing — by the 
following diagram 

A + A ^'^« (1 + e) A (3.1) 

where e > is a constant expressing the increase in species density, and tq is 
a rate constant expressing the velocity of the transformation, at the population 
dynamics time scale. 

Based on the mass action law of §2, the time evolution associated to the 
transformation (3.1) is 

—— = roeAN 

Multiplying by e the second equation in (3.2) and adding to the first one, it follows 
that the time variation of A(t) -|- eA(t) is zero, and, therefore, A(t) + eA{t) = 
constant. With K = A(0) + eA(0), and eliminating A from equations (3.2), we 

obtain 

dN __, _, _/ A~ 



^^ roNiK-N):=rN\l--\ (3.3) 

where K is the carrying capacity, and r = tqK is the intrinsic growth rate of the 
population. The species dynamics (3.3) has the solution A(t) = KN{0)e'^^/{K + 



(e'~* — l)A^(O)), and predictions about the values of the density of a population can 
be obtained by fitting a time series with the explicit solution N{t). In the limit 
t -^ oo, N{t) -^ K. 

In order to make the mechanism (3.1) more realistic, we introduce death rate 
occurring at the ecological time scale: 

A + A^^'^o (l + e)iV 

(3-4) 

N ^'^ ^A 

where the second diagram represents the death of individual with death rate (i, 
and /3 is a recycling constant determined below. The time evolution of the trans- 
formation (3.4) is now 

— - = tqcAN - dN 

dA ^3-') 

-- = -roAN + d(3N 
dt 

Imposing a conservation law of the form, 7A + A^ = constant, and introducing 
(3.5) into the equation 7A + A^ = 0, we obtain, 7 = 6 and /? = 1/e. Therefore, 
with K = N{0) + eA{0), and eliminating A from equations (3.5), we obtain 

dN / A^\ 

-— = roN(K -N)-dN -.^rN il ] - dN (3.6) 

dt \ K J 

where K is the carrying capacity, r = tqK is the intrinsic growth rate of the 
population, and d is the death rate. If r > (i, (3.6) has a stable equilibrium 
solution for N = K{r — d)/r. If r < (i, the only (stable) nonnegative solution 
is A^ = 0. But now, the carrying capacity parameter K is not the equilibrium 
value attained by the population in the limit t — > 00, instead it is the value of 
the conservation law associated to (3.5). In this case, the solution of (3.6) is 
N{t) — e'^^K{d — r)/(Ce'^* — re^*), where C = r + K{d — r)/N{0), and time series 
fitting of observational data is straightforward. 

The recycling condition A^ -^ PA has been introduced in (3.4) in order to 
have a conservation law, which leads to the decoupling of system (3.6), and the 
determination of an explicit solution. But, for example, in bacterial batch cultures, 
where the logistic equation is tested, the recycling condition is not verified. In these 
systems, after the bacterial population has exhausted the resources, the population 
density decreases, a feature that can not be obtained with the recycling condition 
in the logistic equation. 

Dropping the recycling condition in (3.4), we obtain the mechanism 

A + N ^'''-^ (l + e)A^ 

(3-7) 
N ^'^ B 
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where B is some nonrecycling resource. In this case, the evolution equations 
associated to A and A^ are 

—- = tqcAN - dN 

but no conservation law exists enabling its integration. Therefore, the fitting 
parameters are more difficult to estimate. Clearly, the system {A, N} is open, but 
the system {A, N, B} is closed and has a conservation law (eA+N+B = constant) . 

In Fig. 1 we show the graph of the solution N{t) of the logistic type models 
(3.3), (3.6) and (3.8). Comparing the qualitative behavior of the three systems, 
we take the following conclusions. 
i) For large values of available resources A, the solutions of the three systems 

are quantitatively similar, 
ii) The carrying capacity parameter K only coincides with the equilibrium value 

of the population if there are no deaths and populations remain constant after 

exhausting the resources, 
iii) In the exponential growth phase, the three models give qualitatively and quan- 
titatively similar results, 
iv) If the death rate is small compared with the rate constant ro, the maximum 

density of populations calculated by the three logistic models is approximated 

by the carrying capacity K, and the population densities in the exponential 

phase of growth are similar. 

In the following, we will always keep the biomass recycling hypothesis, leading 
to a conservation law and, therefore, to a carrying capacity. 

Specifying intermediate internal states in the life cycle of a species, we now 
show that, under a steady state approximation, we also obtain a logistic equation. 

Suppose that the species A^ has n behavioral internal states, A^i, . . . , A^„, and 
that reproduction occurs according to the mechanism 

A + Ni^N2 Ni^/3iA 



A + N2^Ns N2^ I32A 



A + iV, ^(l + e)iVi Nr^^PnA 



(3.9) 
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where the (3i are constants to be determined later, in order to introduce a conser- 
vation law. The dynamical equations of mechanism (3.9) are 

rlA "^ '^ 

i=l i=l 

^ = -riAiVi + r„(l + e)AN^ - dN, (^'^O) 

at 

—^ = n-iAN^-i - nAN^ -dN,, i = 2, . . . , n 
dt 

We now impose the conservation law, 

dA^ ^^^...^ ^ = (3 11) 

dt dt ^ dt 

Introducing (3.10) into (3.11), and solving this equation for any A and A^^, we 
obtain the parameter values 

'fi = - + i - 1 , Pi ^ 'fi , i ^ 1, . . . ,n (3.12) 

e 

and, from (3.11), the conservation law is 

A + 7iiVi + ... + 7„iV„ = K (3.13) 

We now introduce a steady state assumption, over the internal states, 

dNi 

^ = 0,z = 2,...,n (3.14) 

Solving equations (3.14), we obtain, 

where, by (3.13), A is a function of the N^. But, taking the limit, (i ^ 0, we 
obtain, A^^ = ri-iNi_i/ri, which, by induction, gives, 

A^, = — A^i,z = 2, ...,n (3.16) 

n 

and, with A^ = A^i + . . . + A^^, by (3.10), (3.11), (3.13) and (3.16), 

(3.17) 



which is the logistic equation. 

Therefore, at the ecological scale, internal intermediate states do not introduce 
further dynamical changes in population dynamics equations, as far as the death 
rate d is small. 

For general systems with several basic nutrients, one conservation law for each 
resource should be introduced. 

4. Monod autotrophs 

In order to derive the mechanisms for Monod type population dynamics mod- 
els, we consider a system with an autotroph A^, which can be found in two 
states: searching for nutrient, Ng, and processing (handling) nutrient N^, with 
N = Ns + Nh. When the autotrophs find nutrient, they switch from searching 
to handling, increasing their biomass. At a behavioral time scale, handling au- 
totrophs decay to searching autotrophs, and reproduce. The death rate is the 
same for both handling and searching autotrophs. The kinetic mechanism is thus: 

(4.1) 

where e > is the conversion constant accounting for the increase in species 
density, ri and r2 are the rates at which processes occur, and /? is a recycling 
constant that will be calculated later under the assumptions of a conservation law. 
The resource density is A. 

Let us now apply the formalism of §2 to the interaction given by diagram 
(4.1). By (2.1) and (2.2), we obtain 

dA 

— = -nANs + dP{Ns + Nn) 

^^■^ .^riANs + r^Nh-dN, (4.2) 



dt 

dNh 



ri{l + e)ANs-r2Nh-dNh 



dt 
as evolution equations of resources and individuals. 

From the point of view of population dynamics, we can count species numbers 
but resources are difficult to estimate. Therefore, in order to apply and compare 
the predictions of system (4.2) with a real system, we must be able to rewrite (4.2) 
without modelling explicitly the resource density A. As in the logistic equation, 
the only way to do this is to impose a conservation law, 

dA dNs dNh _ 
dt dt dt 



Introducing (4.3) into (4.2), and solving for the parameters, we obtain, 

7 = e,/? = - (4.4) 

e 

and the conservation law, eA + Ng +Nh = constant = K. Under these conditions, 
the system described by the mechanism (4.1) with j3 given by (4.4) is, 

dN ri 

— ^ = --Ns{K -Ns- Nh) + r2Nh - dN^ 

6 /A r\ 

dt e 

where K = eA + Ns + Nh is the carrying capacity. 

But, at the ecological time scale, our goal is to follow the time evolution of 
the total density of a population, N = Ng + Nh, without knowledge of behavioral 
states of the population. Therefore, with N = Ns + Nh, system (4.5) is rewritten 

as, 

dN 

— = ri{N - Nh){K - N) - dN 

^ = '-^^^±^iN - Nh){K - AT) - r,Nh - dNh 
dt e 

It is natural to assume that there are two time scales in this problem: at the 

ecological time scale the dynamics of Nh is so fast that A^^ is constant. This allows 

us to apply the steady state assumption, 

^ = (4,7) 

dt ^ ' 

implying that, at the ecological scale, 

-^ = nN (K - N) —^, dN := riN(K-N) -^ --dN 

dt l + ^(feS|(^-A^) ^1 + S{K-N) 

(4.8) 
where 6 = ri(l + e)/e{d + r2). System (4.8) has an equilibrium state for A^ — 
{Kri — d{l + 6K))/{ri — dd). This equilibrium state equals the value of the 
carrying capacity K when d -^ 0. 

If we take the behavioral dynamics in (4.1) infinitely fast when compared with 
the ecological time scale, r2 ^ oo, 5 ^ and (4.8) reduces to the logistic equation 
(3.6). For the steady state assumption (4.7) to be valid, the dynamics of Nh must 
be much faster than the dynamics of A^, which is simply obtained for large r2- 

In Fig. 2, we compare density growth curves of the one species model (4.8) 
with the logistic model (3.6), for the same parameters values. The effect of the 
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introduction of an intermediate state, delays the growth and shghtly changes the 
steady state. This contrasts with model (3.9), where the introduction of inter- 
mediate behavioral states corresponding to the different states of the life cycle of 
an individual do not change the overall time behavior of the population at the 
ecological time scale. 

5. Modified Lotka-Volterra trophic chains 

We now consider a kinetic mechanism in accordance with the assumptions 
underlying the Lotka-Volterra prey-predator equations. We represent resources by 
A, prey by Ni and predators by N2. In order to be able to derive a conservation 
law we introduce a recycling mechanism depending upon two unknown parameters, 
f3i and (32- Under these conditions, the prey- predator mechanism is: 



(5.1) 



A + Ni^{l + ei)Ni Ni^/3iA 
Ni + N2^{l + e2)N2 N2 ^ (32A 

where ei and 62 are conversion factors, and di and d2 are death rates. 
The dynamic equations for mechanism (5.1) become, 

dA 

-— = -riANi -F dil3iNi + c?2/92A^2 
dt 

^ = einANi - r2NiN2 - diNi (5.2) 

dt 

dN2 

—— = r2e2NiN2 - d2N2 
dt 

Imposing the conservation law, 

A + 7iA^i + 72A^2 = K (5.3a) 

and after derivation, by (5.2), we obtain 

/^i = 71 = 1/ei , /32 = 72 = 1/(6162) (5.36) 

Using the conservation law (5.3) to eliminate A from (5.2), the equations for 
the prey-predator mechanism become 

"^^^ ein (k-^-^\Ni- r2NiN2 - diNi 



dt \ ei 6162/ /g ^N 

7 TV 7- ^ ' ^ 

—^ = r262A^lA^2 - d2N2 
dt 
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These equations are functionally equivalent to the usual Lotka-Volterra equa- 
tions. However, due to mass conservation, the parameters no longer have exactly 
the same meaning. In fact, prey growth is not only controlled by predation but 
also by the fact that predators retain nutrient within them. 

We can generalize the preceding mechanism to trophic chains of arbitrary 
length: 



ri 



di 1 



A + A^i ^(l + ei)iVi Ni^—A 

ei 

Ni + N2'^{l + e2)N2 



This system has the conservation law: 



N2 


d^ 


' A 

6162 




Nr, 


dr^ 


1 


■A 



(5.5) 



ei...er 



^ + ^ + ^ + ,„ + ^Xl^ = ;, (5.6) 

ei 6162 6162. ..e„ 

Using the conservation law to eliminate A, the dynamical equations become 

—— = eiri A ... Ni - r2iViiV2 - aiiVi 

at \ ei 6162 6162. ..en J /g ^n 

T TV 7- ^ ' ^ 

—-^ = TiCiNi^iNi - diNi , z = 2, . . . , n 
dt 

For trophic chains of length greater than two, the equation for the basal 
species is functionally different from the Lotka-Volterra food chain. The basal 
species is controlled by all the other species, since they are all retaining nutrient 
(this effect increases with the increase in nutrient retained in the nonbasal trophic 
levels). 

6. Trophic chains with internal states 

We now derive the basic population dynamics equations for the time evolution 
of n species in a food chain, assuming that all the species involved have some 
refractory time, during which they are not able to consume resources. 

Consider a food chain with n species, Ni, . . . , N^ and a primary resource 
A. Suppose in addition that each species has two states Nis and Nih, where 
subscripts s and h stand, respectively, for "searching for prey" and "handling for 
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prey" . Introducing this distinction we have A^^ = Nis + Nih . With a characteristic 
time tih, handhng predators finish handhng their prey and return to the searching 
state. Reproduction transforms searching into handhng predactors. Under these 
conditions, the mechanism for the trophic chain is 

A + iVi, ^(l + ei)iVi^ iVis^AA 
,-1 

Nis + N2s ^ (1 + e2)N2h N2h ^ /32A 
Nih + N2s^{l + e2)N2h 

N2H - N2S N2S ^ /32A (6.1) 



Nn-is + Nns ^ (1 + e,,)N^h Nnh ^ (3nA 
N^-ih + A^n. ^ (1 + e^)N^h 

Applying the mass action law of §2, we obtain the evolution equations, 

— = -riANis + 5]] Arf,(Ar,s + N,h) 



^' -- -riAN^s + tih^ih - diNis - r2NisN2s 



dt 

dNih 

dt 



ri(l + ei)ANis - t^h^ih - diNi^ - r2N^hN2s 






dNih 

—— = ri[i. + ei)[i\i-is + i\i-ih)iyis - ^ih 






dN^h 



dt 
where, i = 2,...,n — 1. We now impose a conservation law of the form, 
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where 7^ are unknown parameters. Introducing (6.2) into (6.3), and solving for 7^ 
and (3i, we obtain 

7i = , /?i = 7i , z = 1, . . . , n (6.4) 

ei . . .Ci 

Hence, the conservation law is 

A + —Ni + ... + i A^, = K (6.5) 

d ei...en 

where K is the carrying capacity, and Ni = Ni^ + Nih . 

Let us now introduce a steady state assumption over the behavioral states 

N^h■. 

-_ =0,2 = l,...,n (6.6) 

With, Ni = Nis + Nih, and writing the system of equations (6.2) as a function of 
Nn and A^ns, we obtain the dynamic equations for species densities, 

^ = eiri(K - -N^ - . . . ^ N^)N^, - r^N^N^s - diN^ 

at ei ei . . . e^ 



—^ = e,r^N,sN,-i - r^+iN,N,+i, - d,N, (6.7) 

at 



at 

In order to eliminate Nis from (6.7), we solve the system of equations (6.6) together 
with the relation A^^ = Ni^ + Nih. Therefore, (6.7) together with (6.6) are the 
general equations for a trophic chain. 

Let us now analyze the case n = 2. In this case, (6.7) and (6.6) simplify to 

^ = emiK - -Ni - J-N2)N^s - r2mN2s - diN^ 
at ei 6162 

dN2 



dt 

dNi 



h 



e2r2NiN2s - d2N2 

^^ nil + ei)AN,s - t-^{N^ - Nis) - di{Ni - N^^) - r2{m - ms)N2s = 

^ = r2(l + e2)AriAr2, - t-^{N2 - N2s) - rf2(A^2 - N2s) = 

(6.8) 
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Solving the two last equations in (6.8) for A^i^ and A^2s, we obtain, 



Nos = N, "^^ 



W+^2+r2(l + e2)iVi 



tlh +di+ri{l + ei)(A' - j-Ni - ji-iVj) + r^N-, -,'^*^;^ ,„ 

->-"■ ^ / \ ei 6162 ' t^^ ^c(2+r2(l+e2)A'i 



Therefore, system (6.8) simplifies to 

dNi 

~dr 

dN2 
dt 



ein{K 


1 

ei 


iVi-- 


1 

2162 


e2r2N2s 


iVi- 


-d2N2 





iV2)iVis - r2NiN2s - diN, 



(6.9) 



(6.10) 



where Nis and A^2s are given by (6.9). 

Equation (6.10) describes a trophic chain with a prey, a predator and a re- 
newable resource. The carrying capacity of the system is the constant K. The 
introduction of the intermediate states Nih, specifying the existence of a refrac- 
tory time where species are not consuming resources, implies a functional response 
of HoUing type II, with ecologically meaningful parameters. If we consider that 
the basal species in A^i has negligible handling time, t^^ -^ 00, then, by (6.9), 
Nis -^ A^i, and (6.10) reduces to 

^ = eiri{K - -Ni - -^N2)N, - r2NiN2—^ /'^ ^,f' - diNi 

dt ei 6162 t^h +d2 + r2il + e2)Ni 



e2r2NiN2—, , ^^ ' 7^ -— - rf2A^2 



1 = e2r2NiN2^. — = 

dt t^^ +d2+r2il + e2)Ni 

(6.11) 
wnich is a modified Rosenzweig-MacArthur (1963) predation model. 

In the case the handling times t^^ and t^^ go to infinity, by (6.9), Nis con- 
verges to Ni, and (6.10) reduces to the modified Lotka-Volterra equation (5.4). 

We now compare the time evolution of system (6.10) with the modified prey- 
predator Lotka-Volterra system (5.4). We take the parameter values: ri = 1.0, 
r2 = 0.2, ei = 1, 62 = 1, di = 0.1, ^2 = 0.2, t'^ = 1.0, t'^ = 2.0, K = 8, 
Ni{0) = 1.0 and A^2(0) = 0.1. In Fig. 3, we depict the time evolution of the two 
models, for prey and predators. In the case of model equations (5.4), the system 
attains an equilibrium value given by A''! = d2/e2r2 and N2 = {r2{Keiri — di) — 
d2Ti)/{r2{ri + 62^2)). However, for system (6.10), we obtain stable oscillations 
corresponding to a limit cycle in phase space. 
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For the general case of an arbitrary trophic chain we have equation (6.2) 
together with the carrying capacity relation (6.5), and, eventually, the steady 
state conditions (6.6). 

7. Conclusions 

We have developed a systematic formalism, based on chemical kinetics, for 
the derivation of equations in population dynamics based on the mechanisms of 
interaction between individuals. We have started from the simplest equation, the 
logistic, and then introduced successive levels in a trophic chain. The dynamical 
equations are derived using the laws of mass action and mass conservation, and, 
when necessary, a steady state assumption. Our approach includes any kind of 
trophic interaction between species and resources and internal states. This ap- 
proach gives a consistent mechanistic basis for the derivation of the trophic chain 
equations of population biology, making it possible to settle several controversies in 
Ecology. Moreover, this formalism allows the precise development of more compli- 
cated models, with the introduction of more mechanisms and interactions, allowing 
the development of extensions of the logistic equation with precise applicability 
conditions. 

This approach has the advantage that all the parameters have an a priori 
biological meaning. Moreover, we show the necessity of the introduction of a con- 
servation law relating populations and primary resource densities. In this context, 
the carrying capacity parameter is the value assumed by the conservation law, 
and appears in all population dynamics models, from one species to arbitrary food 
webs. However, the numerical equilibrium value attained by the populations only 
equals the carrying capacity in the case of the logistic equation. 

The use of the mass action law, together with its microscopic foundation, 
allows the development of models with stochastic fluctuation around a mean value 
(Bailey, 1964; Haken, 1983), as well as the development of equations with different 
growth forms and the introduction of individual variability. 

Finally, one of the consequences of the systematic use of the mass action law 
is the proof of the empirical observation that the growth rate r, and the carrying 
capacity K in the logistic equation, both increase proportionally to enrichment. 
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Figure Captions 
Figure 1: Solution N{t) of the logistic type models (3.3), (3.6) and (3.8). The 
parameter values are: e = 1, tq = 1, (i = 0.2, K = 20, A^(0) = 0.1 and A{0) = 
{K-N{0))/e = 19.9. 

Figure 2: Solution N{t) for the logistic and Monod type model (4.8). The pa- 
rameter values are: e = 1, ri = 1, r2 = 10, d = 0.2, K = 20, A^(0) = 0.1, 
A{0) = {K - A^(0))/e = 19.9 and 5 = 0.196. In the limit 5 ^ (or r2 ^ oo) the 
solution of the Monod type model approches the solution of the logistic equation. 

Figure 3: Solutions N{t) of the modified Lotka-Volterra model (5.4) and trophic 
model (6.10). The parameter values are: ri = 1.0, r2 = 0.2, ei = 1, 62 = 1, 
rfi = 0.1, d2 = 0.2, q^ = 1.0, t^^^ = 2.0, K = 8, iVi(O) = 1.0 and A^2(0) = 0.1. 
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This figure "fig 1. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/adap-org/9902002vl 



This figure "fig2.GIF" is available in "GIF" format from: 
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